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ABSTRACT:  It  is  shown  theoretically  that  the  pressure 

rise  across  the  shock  wave  about  a body  of  revolution 
decays  asymptotically  as  the  inverse  three -fourths  power 
of  the  radial  distance  from  the  body  axis.  Experimental 
data  is  presented  to  support  the  theory. 
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ON  TEE  ATTENUATION  OF  THE  SHOCK  NATE 
ABOUT  AN  AXIALLT  - SYMMETRIC  BOOT 


INTRODUCTION 

1.  The  attenaution  of  the  shock  wave  about  a two-dimensional  body 
Is  discussed  by  Munk  and  Crown  in  reference  (a).  It  is  shown  that  the 
difference  between  the  shock  angle  and  the  free -stream  Mach  angle  decays 
inversely  as  the  square  root  of  the  distance  away  from  the  body  axis. 
Hence,  the  shock  wave  contour  at  large  distances  from  the  body  becomes 
parabolic,  having  no  asymptote.  In  addition  to  practical  considerations 
such  as  the  extrapolation  of  a given  portion  of  a shock  wave  to  obtain 
wave  drag,  as  Indicated  in  reference  (a),  the  question  of  whether  or  not 
a shock  wave  about  an  axially- symmetric  body  has  an  asymptote  and  the 
nature  of  the  shock  decay  are  of  fundamental  importance . 

2.  In  19^9  the  author  was  interested  in  the  problem  of  the 
antentuation  of  shock  wave  about  a body  of  revolution  and  in  discussion 
with  various  people  was  verbally  Informed  that  some  work  had  been  done 
on  this  problem  which  showed  that: 


* where  x-  streamwlse  coordinate 

r = radial  distance  from  body  axis 
Mi  = free -stream  Mach  number 
a = an  arbitrary  constant 

3.  The  work  containing  the  formula  as  well  as  the  derivation  could 
not  be  located  in  the  literature.  A preliminary  check  of  this  formula 
with  the  then  available  experimental  data  did  not  yield  satisfactory 
agreement.  In  order  to  evaluate  the  correctness  of  the  formula  the 
author  of  this  report  made  an  independent  analysis  which  is  presented 
herein.  The  results  yielded  a formula  essentially  the  same  as  given 
above.  Substantiation  of  the  formula  was  obtained  when  recent  accurate 
and  reliable  data  were  used. 

h.  Thi3  report  contains  the  author’s  derivation  of  the  theoretical 
formula  for  the  attenuation  of  the  shock  wave  about  a body  of  revolution 
(or  any  three-dimensional  body,  for  that  matter)  and  a presentation  of 
experimental  data  to  substantiate  it. 


ANALYSIS 

The  Characteristic  Equations: 

5.  If  the  shock  wave  curvature  is  small,  as  it  is  as  large  distances 
from  a body,  then  with  little  error  we  may  take  the  flow  behind  the  shock 
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to  be  lsentroplc.  The  characteristic  equations  then  take  the  form 


dd  + dB  7 


C*n  0 c«v»  ol 

£<  n (pit  d 


djf  XT  Q 


(2) 


where 


0 is  the  expansion  angle 

0 is  the  flow  inclination 

01  is  the  Mach  angle 

V*  is  the  radial  distance  from  the  body  axis , and 
the  upper  and  lower  signs  refer  respectively  to  families  I and  II  charac- 
teristics as  indicated  in  Figure  1.  When  the  difference  from  free-stream 
conditions  is  small , the  characteristic  equations  may  be  approximated  by 


LO7/LO  7 O ^ = 0 

Z 

6.  From  the  geometry  of  Figure  1,  It  can  be  seen  that 

(at ~d)  Zt**  ~ O') 

A £ (ol  Id  ) &V,  (d  t<p * d ) 

~ (*-  *t)  t Q - € 

’t-'i  C.V,  2 <u, 


where  £ ~ ^ 


(5) 


C is  the  shock  wave  angle 


and  the  subscript  "1"  refers  to  free-stream  conditions. 


(6) 


(7) 


and  from  the  condition  of  is  entropy  lasted  lately  behind  the  shock  wave 
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ve  have  the  relation 


(8) 


From  the  veak  oblique  shock  formulas  ve  find 


(9) 


Substituting  thesr  formulas  into  equation  (5)>  ve  obtain 


^ J L*  let, 

Thus  it  can  be  seen  that  is  one  order  of  magnitude 

Ignoring  in  the  characteristic  equations  for  the 


(10) 


smaller  than  . 
II  f amily,  ve  find 


/ 


(11) 

or 

(8) 

at  any  point  behind  the  shock  vhere  conditions  are  not  far  from  free-Btream 
values.  Substituting  equation  (8)  into  the  characteristic  equation  for 
family  I,  ve  obtain 


1 (JlQ  y 


or  upon  integration 

0 r 
&' 


0 d/ - 0 

? 


(12) 


(13) 


vhere  the  primed  values  represent  the  constants  of  integration.  This  equa- 
tion is  the  one  upon  vhich  the  asymptotic  nature  of  the  shock  vave  is  based. 

The  Shock  Wave: 

7.  The  slope  of  the  Mach  lines  of  family  I is  given  by 

T -/a,-,  d (it) 

h 
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where 


4>  r <*V  Q 


For  small  deviations  from  free-stream, 


we  can  write,  as  before, 


(15) 


(16) 


where  we  now  take  X,  / to  be  the  origin  of  the  expansion  waves  which 
weaken  the  shock.  This  point  may  be  real  or  imaginary . For  a cone-cylinder 
model,  this  point  is  at  the  shoulder.  However,  the  actual  existence  of  such 
a point  is  immaterial  to  asymptotic  behavior  of  the  shock  wave. 

8.  Using  the  first  two  terms  of  a Taylor's  lxpansion,  we  may  write 
the  slope,  ( 14) , as 


; v " ^ % r-~  j 

r 4moIi  / / 4 JL  (i/>  -tf,  ) j - / 

7%  i J K J 


(17) 


Let 


(18) 


which  is  constant  along  a given  Mach  line.  Then  equation  (17)  becomes 
(Xy*  - 4ou^dt,  cW  - ^ JLf  (19) 


K 4 

t: 

K 


and  upon  integration,  yields  the  equation  of  the  shock  wave: 


K+r? 


- df 


(20) 


Where  the  subscript  "s"  refers  to  conditions  immediately  behind  the  shock 
wave.  Nov  the  slope  of  the  shock  wave  is  given  by 

z 4am(oI,+£)  (21) 

Jy, 
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r cJ-4,  -Alfa 


Differentiation  of  aquation  (20)  yields  another  expression  for  the  slope 
•f  the  shock  ware:  . , . , i P % 


Combining  with  (22)  we  have 

g - /?(  ^ AX 


9 •> 

7^ 


* 


in  I 


M'-i 


Fro*  (18) 


5-  — Ms 
dx  \ K*Ss  J 


K.  = ^ /£' 

/*<£/ 

or  v, 

K*  4.  (*s -<M 

Mi/ 

From  the  foregoing  equations,  it  can  he  shown  that 


</s  - <S,  'ZT'S 


and  therefore 


k-  2Jt‘  efs{ 

Mm 


Thus  the  quantity 

K_ 


/M^t 
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Differentiating  equation  (27)  ve  obtain 

djC  __  l M,  ^ ^cf  | 


(29) 


The  Integral  in  equation  (24)  any  be  evaluated  an  follows: 

Jr  *<«W  J,‘  (k<S?)v 

_ l [(&/*■)  ■ XK.  jL>  (Xi  ft).  //  / 2r  £ j1 

' *vi/?  fti  V 


Substituting  these  relations  into  aquation  (24)  and  dropping  the  subscript 
"s\  ve  obtain  for  the  equation  of  the  shock  vave  the  relation 

V f l r 

3/m^i  r* 


A r-  i ,£ 


/-  i ^ f X r 

/5<T/  /' 


(31) 


This  differential  equation  can  be  satisfied  in  the  region 

only  with  £ lu  f > Q 

M,1 

The  quantity  remains  finite  for  / ( 4^  ( o0  . 

Considering  f as  a higher  order  infinitesimal  and  Integrating 

/ft,  f* 

equation  (31) , ve  obtain  a formula  for  the  asymptotic  shape  of  the  shock 
vave:  , 

£*  J- 


rv¥  <*> 

vhere  k is  the  constant  of  integration  vhlch  may  be  taken  equal  to 

fa/W 

and  X))  is  a point  at  vhlch  S Is  known. 
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9.  Equation  (32)  is  the  one  ve  have  been  seeking.  It  can  be  seen  to 
be  equivalent  to  equation  (1)  by  differentiating  the  latter. 

10.  From  equation  (1)  or  (32) , it  follovs  that  vhile  the  shock  angle 
approaches  the  free-stream  Mach  angle  as  f-y-co  , the  shock  wave  contour 
itself  does  not  have  an  asymptote. 


11.  Furthermore,  from  equations  4.27  and  4.29  of  reference  (b),  it 
follovs  that  the  fractional  pressure  rise  across  a veak  shock  wave  is 
given  by 


4k 

p 


(33) 


This  pressure  disturbance  decays  also  as  the  inverse  three-fourths  power 
of  the  distance  from  the  body  axis. 

12.  In  order  to  demonstrate  the  validity  of  the  three-fourths  power 
lav,  equation  (32),  it  vould  seem  logical  to  plot  experimental  values  of 
f on  log- log  graph  prfper.  The  theoretical  relation,  which  is  the 
asymptotic  solution  for  f-yoo , then  appears  as  a straight  line  with  a 
negative  slope  of  3/4.  This  comparison  is  given  in  the  following  section. 

COMPARISON  WITH  EXPERIMENT 

13.  The  data  presented  herein  were  obtained  from  spark  shadowgraphs  of 
firings  made  in  the  NOL  Pressurized  Ballistics  Range.  The  velocity  was 
determined  accurately  electronically  and  the  sound  speed  obtained  from 
temperature  measurements;  thus  the  Mach  number  could  be  calculated  accu- 
rately. It  is  believed  that  the  greatest  source  of  error  was  incurred  in 
the  measurement  of  the  shock  angles  which  was  done  on  Bausch  and  Lomb 
comparators.  The  scatter  of  the  data  is  thought  to  be  indicative  of  the 
measuring  accuracy  and  this  appears  to  be  within  satisfactory  limits. 

14.  Figure  2 shows  a comparison  of  theory  and  experiment  for  the  shock 
wave  about  a sphere.  The  agreement  can  be  seen  to  be  very  good.  Figure  3 
shows  a similar  comparison  at  a higher  Mach  number.  It  should  be  emphasized 
that  equation  (32)  is  an  asymptotic  solution  and  hence  the  divergence  of  the 
theoretical  curve  from  the  experimental  data  at  small  distances  from  the  body 
is  not  contradictory.  At  large  distances  from  the  body  the  asymptotic  solu- 
tion appears  to  be  approached  very  closely. 

15.  Figure  4 illustrates  the  attenuation  of  the  shock  wave  for  a missile. 
The  theoretical  curve  again  agrees  very  well  with  the  experimental  data.  Some 
of  the  irregularities  at  large  radii  may  be  due  to  the  influence  of  the  fins. 

16.  The  data  for  another  missile  is  presented  in  Figure  5*  Here,  too, 
the  agreement  with  theory  is  close  up  to  about  15  body  radii  away  from  the 
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missile  axis  at  which  point  the  shock  strength  appears  to  remain  constant. 
This  effect  Is  believed  due  to  yaw  of  the  missile.  That  this  is  the  expla- 
nation might  have  been  demonstrated  by  measuring  the  shock  angles  for  the 
other  side  of  the  missile.  However,  in  order  to  obtain  shock  waves  which 
extend  to  as  large  a distance  from  the  body  as  possible,  plates  were 
selected  on  which  the  missile  was  off-center.  Hence  no  such  comparison 
could  be  made.  It  seems  obvious,  however,  that  this  constant  shock  angle 
cannot  be  maintained  indefinitely  if  the  drag  is  to  remain  finite.  Thus, 
when  the  expansion  waves  from  the  tail  of  the  missile  begin  to  reach  the 
shock,  it  will  begin  to  decay  once  more  as  the  inverse  three-fourths  power 
of  the  distance. 


CONCLUSION 

17.  It  is  believed  that  the  experimental  data  presented  herein  confirm 
the  basic  validity  of  equation  (32)  for  the  asymptotic  attenuation  of  the 
shock  wave  about  a body  of  revolution. 


FURTHH?  HISTORICAL  NOTE 

18.  After  the  completion  of  this  report  and  just  prior  to  publication 
several  pertinent  references  were  brought  to  the  attention  of  the  author. 

•Hie  unspecified  reference  alluded  to  in  the  introduction  was  found  to  be 
the  work  of  Whitham,  reference  (c).  Therein,  a derivation  of  the  formula 

in  question  is  given.  However,  prior  to  him,  DuMond  et  al  in  reference  (d) 
gave  a different  derivation  and  presented  experimental  data  for  large 
distances  from  the  missile.  Furthermore,  as  far  as  the  author  has  been 
able  to  determine,  the  formula  was  first  mentioned  by  Landau  in  reference  (e) 
where  only  a simple  statement  of  the  formula  is  given. 

19.  The  present  report  with  new  experimental  data  was  considered  to 
be  of  sufficient  interest  to  justify  publication  in  addition  to  the 
preceeding  reports. 
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